RF signal processing in multi-antenna systems

ABSTRACT

A method for antenna subset selection by joint processing in RF and baseband in a multi-antenna systems. L t  input data streams are generated in a transmitter for either diversity transmission or multiplexing transmission. These streams are modulated to RF signals. These signals are switched to the t branches associated with the t transmit antennas, and a phase-shift transformation is applied to the RF signals by a t×t matrix multiplication operator Φ 1 , whose output are t≧L t  RF signals. These signals are transmitted over a channel by t antennas. The transmitted signals are received by r antennas in a receiver. A phase-shift transformation is applied to the r RF signals by a r×r matrix multiplication operator Φ 2 . L r  branches of these phase shifted streams are demodulated and further processed in baseband to recover the input data streams.

FIELD OF THE INVENTION

The invention relates generally to multi-antenna communication systems, and more particularly to processing RF signals transmitted through spatially correlated channels in such systems.

BACKGROUND OF THE INVENTION

Multi-antenna systems are widely considered to be the most viable way for significantly increasing the bandwidth efficiency of wireless data transmission systems. In MIMO (multiple input multiple output) systems, multiple antennas are deployed both at the transmitter and the receiver. In MISO (multiple input single output) systems, the receiver has only one antenna, and the multiple transmit antennas are used for transmit diversity. In SIMO (single input multiple output) systems, the transmitter has a single antenna, and multiple antennas are used at the receiver.

Given multiple antennas, the spatial dimension of the channel can be exploited to improve the performance of the wireless link. The performance is often measured as the average bit rate (bit/s) the wireless link can provide, or as the average bit error rate (BER), depending on the application.

Given a multi-antenna channel, a duplex method, and a transmission bandwidth, the multiple-antenna system can be categorized as narrowband or wideband (i.e., the channel is flat or frequency selective fading within the system bandwidth), and possessing either full, partial or no channel state information (CSI).

Multiple antennas provide enhanced performance. Studies show that the multiple antennas can be used to provide spatial diversity, and/or can increase the information-theoretic capacity (data rate), see, e.g., Vaughn and Anderson, Channels, propagation, and antennas for mobile communications, IEE Press, 2003, Telatar et al., “Capacity of Multi-Antenna Gaussian Channels,” European Trans. on Telecomm., Vol. 10, No. 6, pp. 585-596, November-December 1999, Winters, “On the Capacity of Radio Communication Systems with Diversity in Rayleigh Fading Environments,” IEEE J. Selected Areas Comm., 1987, and Tarokh et al., “Space-Time Codes for High Data Rate Wireless Communication: Performance Criterion and Code Construction,” IEEE Trans. Inform. Theory, Vol. 44, pp. 744-765, March 1999.

However, as a disadvantage, operation with multiple RF chains increases complexity and cost. An antenna selection technique can be used to determine an optimal subset of antennas. This still yields improved performance while reducing the number of required RF chains, see Molisch et al., “Capacity of MIMO Systems with Antenna Selection,” Proc. IEEE Intl. Comm. Conf., pp. 570-574, 2001, and Gore et al., “MIMO Antenna Subset Selection with Space-Time Coding,” IEEE Trans. Signal Processing, Vol. 50, No. 10, pp. 2580-2588, October 2002.

Most prior art antenna selection techniques have focused only on selecting a subset of antennas before down-conversion and processing in baseband. This (spatial) antenna selection techniques work reasonably well under some circumstances, namely in the case that (i) the number of selected antennas is only slightly smaller than the number of available elements, and (ii) the MIMO channel is spatially uncorrelated.

However, in reality, correlated scattering at both the transmitter and receiver antenna arrays is more usual. Due to the directional transmission in the wireless environment, the signal waveforms can be highly correlated depending on the departure and arriving angles at the antennas.

The prior art antenna selection techniques that operate only in the spatial domain can have a significant degradation in performance when dealing with this correlation. Furthermore, even for weak correlations, the prior art antennas also show a significant performance loss when the number of selected antennas is considerably smaller than the number of available antennas.

Therefore, it is desired to provide an RF signal processing technique that overcomes the problems of the prior art.

SUMMARY OF THE INVENTION

The invention provides a system and method for antenna subset selection in a multi-antenna communication system based on the joint signal processing in radio frequency and baseband. The system is designed for either full diversity transmission or multiplexing transmission, in both correlated and uncorrelated channels.

The invention uses phase-shift operations in radio-frequency (RF) chains. The operations can be performed after antenna selection in the transmitter, or before antenna selection in the receiver, or both.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a block diagram of a multi-antenna system that uses joint RF/baseband signal processing for antenna selection according to the invention; and

FIGS. 2 and 3 are graphs comparing RF signal processing techniques.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT

FIG. 1 shows a multi-antenna system 100 that uses the antenna selection method according to the invention. In a transmitter 1, L_(t) data streams 101 are generated. These streams can be either different signals obtained through a space-time coder, i.e., in multiplexing transmission, or the same signal encoded with different weights, i.e., in diversity transmission. The streams are modulated 110 to passband RF signals. A selection switch 120 links these signals to L_(t) of the t branches 121 associated with the t transmit antennas (t≧L_(t)) 131. In the new invention, the t passband branches 121 are transformed 130 by a t×t matrix operation Φ₁ before they are applied to the t transmit antennas, and passed through the channel 140.

At a receiver 2, the signals are received via the channel 140 by r receive antennas 151. The received signals are transformed 150 by an r×r matrix operation Φ₂, and L_(r) of the transformed signals are selected 160, where L_(r)≦r. The selected signals 161 are demodulated 170 and further processed 180 in baseband for detection of the data streams 101. The concatenation of the L_(t) out of t switch and the t×t matrix operation Φ₁ can be viewed, and implemented, as a t×L_(t) matrix operation in both the transmitter and the receiver.

In order to reduce the implementation cost in radio frequency, the entries of the matrix operations usually are constrained to be pure phase shifters without power amplifiers. They can also be done just in the transmitter, or just in the receiver, e.g., by making the corresponding matrix at the other link end, Φ₂ or Φ₁, identity matrices.

System and Channel Model

For the description of the operating principles of the invention, we use a well known channel model, see Shiu et al., “Fading Correlation and Its Effect on the Capacity of Multi-element Antenna Systems,” IEEE Trans. Commun., Col. 48, pp. 502-513, March 2000, Bölcskei et al., “Performance Analysis of Space-Time Codes in Correlated Rayleigh Fading Environments,” Proc. Asilomar Conf. Signals, Syst. Comput., pp. 687-693, November 2000, Gore et al., “MIMO Antenna Subset Selection with Space-Time Coding,” IEEE Trans. Signal Processing, Vol. 50, No. 10, pp. 2580-2588, October 2002. However, we note that the invention is in no way dependent on the channel model.

The model for the channel transfer function H is expressed as: H=R ^(1/2) WT ^(1/2),  (1) where W is a Rayleigh fading matrix with i.i.d. complex Gaussian entries ˜N_(C)(0,1). The correlation matrices R and T denote receive and transmit correlations, respectively. The sizes of R and T are r×r and t×t, respectively. This model is usually valid when assuming independent transmit and receive correlations.

To facilitate the performance analysis, we introduce the singular value decomposition (SVD) of the transfer function matrix H: H=UΣV*, where U and V are unitary matrices representing the left and right singular vector spaces of H, respectively; and Σ is the diagonal matrix consisting of all the singular values of H. For convenience we denote λ_(i)(H) as the i-th largest singular value of matrix H, and {right arrow over (u)}_(i)(H), {right arrow over (v)}_(i)(H) are the left and right singular vectors of H with respect to λ_(i)(H). Here we use ‘*’ to denote the conjugate transpose of a vector or matrix (It reduces to the complex conjugate for a scalar).

Diversity Transmission

Full-Complexity (FC) MRT/MRC Scheme

The system 100 with diversity transmission can be expressed by a linear equation: {right arrow over (x)}(k)=H{right arrow over (v)}s(k)+{right arrow over (n)}(k)  (3) where s(k)∈C is the transmitting stream, {right arrow over (x)}(k)∈C^(r) is the set of sample stacks of the complex-valued receiver data sequence. The total transmission power is constrained to P. The thermal noise {right arrow over (n)}(k)∈C^(r) is a white i.i.d Gaussian random process with independent real and imaginary parts and variance σ_(n) ²I_(r), and {right arrow over (v)} is a t-dimensional transmitter weighting vector satisfying ∥{right arrow over (v)}∥=1. At the receiver, the received signals {right arrow over (x)} are weighted (W) 152 with complex weights {right arrow over (u)}* and summed (Σ) 153, to give a soft estimate of the transmitted symbol stream.

Without antenna selection, the information stream s(k) is estimated by performing a linear combination of all the r observation streams with a coefficient vector {right arrow over (u)}* at the receiver by: ŝ(k)={right arrow over (u)}*H{right arrow over (v)}s(k)+{right arrow over (u)}*{right arrow over (n)}(k).  (4) The estimate SNR after combining is

$\begin{matrix} {{{a.\mspace{14mu}\frac{{ɛ\left\lbrack \left. {{\overset{\rightarrow}{u}}^{*}H\;\overset{\rightarrow}{v}{s(k)}} \right\rbrack  \right.}^{2}}{{{{\overset{\rightarrow}{u}}^{*}{\overset{\rightarrow}{n}(k)}}}^{2}}} = {\rho\frac{{{{\overset{\rightarrow}{u}}^{*}H\;\overset{\rightarrow}{v}}}^{2}}{{{\overset{\rightarrow}{u}}^{*}}^{2}}}},} & (5) \end{matrix}$ where the nominal SNR is ρ=P/σ_(n) ². To maximize the estimate SNR, maximal ratio transmission and maximal ratio combining are used, i.e., {right arrow over (u)}({right arrow over (v)}) are the singular vectors in U(V) corresponding to the largest singular value λ₁(H). The resulting estimated SNR is then ρλ₁ ²(H).

Conventional Hybrid Selection (HS) MRT/MRC Method

When L_(r) out of the r receive antennas are selected and combined, each selection option corresponds to a L_(r)×r selection matrix on the transfer function, which extracts L_(r) out of the r rows of the matrix H that are associated with the selected antennas. We denote S_(r) as the set of all such selection matrices. Similarly, an L_(t) out of t selection in the transmit end can be represented by a t×L_(t) matrix, and the set of the transmit selection matrices are denoted by S. For any selection option (S₁∈S_(t),S₂∈S_(r)), an optimal SNR is achieved via a similar MRT and MRC on the L_(t) transmit and L_(r) receive branches respectively:

$\begin{matrix} {{\max\limits_{\overset{\rightarrow}{u} \in C^{L_{r}}}{\max\limits_{{\overset{\rightarrow}{v}} = 1}{\rho\frac{{{{\overset{\rightarrow}{u}}^{*}S_{2}{HS}_{1}\;\overset{\rightarrow}{v}}}^{2}}{{\overset{\rightarrow}{u}}^{2}}}}} = {{{\rho\lambda}_{1}^{2}\left( {S_{2}{HS}_{1}} \right)}.}} & (7) \end{matrix}$ The HS-MRC scheme selects the optimal antenna subset selection matrices (S₁, S₂) among all the elements in S_(t) and S_(r) that maximize the estimate SNR

$\begin{matrix} {{{above}\text{:}\mspace{14mu}{SNR}_{HS}} = {\max\limits_{S_{2} \in S_{r}}{\max\limits_{S_{1} \in S_{t}}{{{\rho\lambda}_{1}^{2}\left( {S_{2}{HS}_{1}} \right)}.}}}} & (8) \end{matrix}$

FFT-Based Method (FFTS)

To cope with the highly correlated MIMO channels, the operation (Φ_(t) or Φ_(r)) 150 can be a fast Fourier transform matrix of the form:

${{i.\mspace{14mu}\Phi_{r}} = {\frac{1}{\sqrt{r}}\begin{bmatrix} 1 & 1 & \ldots & 1 \\ 1 & {\mathbb{e}}^{{- j}\;\omega_{r}} & \ldots & {\mathbb{e}}^{{- {j{({r - 1})}}}\omega_{r}} \\ \vdots & \vdots & ⋰ & \vdots \\ 1 & {\mathbb{e}}^{{- {j{({r - 1})}}}\omega_{r}} & \ldots & {\mathbb{e}}^{{- {j{({r - 1})}}^{2}}\omega_{r}} \end{bmatrix}}},$ where W_(r)=2π/r. The matrix is normalized to preserve the noise level. Via a similar MRT/MRC on the selected branches after the FFT, the optimal estimate SNR is:

$\begin{matrix} {{SNR}_{FFTS} = {\max\limits_{S_{2} \in S_{r}}{\max\limits_{S_{1} \in S_{t}}{{{\rho\lambda}_{1}^{2}\left( {S_{2}\Phi_{r}{H\Phi}_{t}S_{1}} \right)}.}}}} & (9) \end{matrix}$

Phase-Shift & Selection Method (PSS)

A more general design in our invention is to allow phase-shift-only entries in the transform Φ₁ 130 (and/or Φ₂ 150). In this case, the transform matrix 130 (150) and the selection 120 (160) can be integrated into one t×L_(t) (or L_(r)×r) matrix with phase-shift-only entries. The set of all such transmit (or receive) matrices is denoted by F_(t) (or F_(r)). In this method, with the two phase-shift matrices F₁∈F_(t) and F₂∈F_(r), the linear combining is performed on the virtual channel F₂HF₁. As above, the SNR for the optimal choice of (F₁, F₂) is

$\begin{matrix} {{SNR}_{PSS} = {\max\limits_{\overset{\rightarrow}{u} \in C^{L_{r}}}{\max\limits_{\overset{\rightarrow}{v} \in C^{L_{t}}}{\max\limits_{F_{1} \in F_{t}}{\max\limits_{F_{2} \in F_{r}}{\rho{\frac{{{{\overset{\rightarrow}{u}}^{*}F_{2}H\; F_{1}\;\overset{\rightarrow}{v}}}^{2}}{{{{\overset{\rightarrow}{u}}^{*}F_{2}}}^{2}{{F_{1}\overset{\rightarrow}{v}}}^{2}}.}}}}}}} & (11) \end{matrix}$ When L_(r)≧2 and L_(t)≧2, the optimal solution for equation (10) can be derived in a closed form. Denote the elements in F₂ as [F₂]_(m,n)=exp(jφ_(m,n)) and the singular vector of H as {right arrow over (u)}₁*(H)=[β₁₁e^(−jφ) ¹¹ . . . β_(r1)e^(−jφ) ^(r1) ], an optimum choice of the phases in F₂ is

$\begin{matrix} \begin{matrix} {{{i.\mspace{14mu}\phi_{1,i}} = {{- \varphi_{i,1}} - {\cos^{- 1}\frac{\beta_{i,1}^{2} + {\beta_{\max}\beta_{\min}}}{\beta_{i,1}\left( {\beta_{\max} + \beta_{\min}} \right)}}}};} \\ {\mspace{31mu}{\phi_{2,i} = {{- \varphi_{i,1}} + {\cos^{- 1}{\frac{\beta_{i,1}^{2} - {\beta_{\max}\beta_{\min}}}{\beta_{i,1}\left( {\beta_{\max} - \beta_{\min}} \right)}.}}}}} \end{matrix} & (11) \end{matrix}$ where

$\beta_{\max} = {{\max\limits_{1 \leq i \leq r}{\beta_{i1}\mspace{14mu}{and}\mspace{14mu}\beta_{\min}}} = {\min\limits_{1 \leq i \leq r}{\beta_{i1}.}}}$ The other rows of F₂ can be arbitrarily designed. The optimum phase shifters in F₁ are derived in a similar fashion with {right arrow over (u)}₁*(H) replaced by {right arrow over (v)}₁*(H). With this choice, the PSS can achieve the same SNR as the FC-MRT/MRC.

When either L_(r)=1 or L_(t)=1 (or both), a closed-form solution to (10) in generally does not exist; we provide a sub-optimal solution which closely approximates the maximization in (10): φ_(1i)=−φ_(r1).  (12)

Multiplexing Transmission

Full-Complexity (FC) Scheme

The system 100 with multiplexing transmission is expressed by {right arrow over (x)}(k)=H{right arrow over (s)}(k)+{right arrow over (n)}(k),  (13) where s(k)∈C^(r) is now a t-dimensional vector denoting the different transmitting sequences. The assumptions on the channel and noise are the same as the diversity transmission case. To support multiple data streams, capacity is an important measure for the capacity of the system.

Without antenna selection, the information rate supported by the original channel can be maximized using a “water-filling” power allocation. The capacity of the original system is

$\begin{matrix} {{C_{FC} = {\sum\limits_{i = 1}^{\min{({t,r})}}\;\left\lbrack {\log_{2}\frac{\rho}{t}{{\mu\lambda}_{i}^{2}(H)}} \right\rbrack^{+}}},} & (14) \end{matrix}$ where [a]⁺ is defined as max(a,0 and μ is the constant satisfying

${\sum\limits_{i = 1}^{\min{({t,r})}}\left\lbrack {\mu - \frac{t}{{\rho\lambda}_{i}^{2}(H)}} \right\rbrack^{+}} = {t.}$

Conventional Hybrid Selection (HS) Method

When L_(t) (L_(r)) out of the t (r) transmit (receive) antennas are selected, the maximum capacity offered by the conventional antenna selection scheme is

${C_{HS} = {\max\limits_{S_{1} \in S_{t}}{\max\limits_{S_{2} \in S_{r}}{\sum\limits_{i = 1}^{\kappa{({S_{2}{HS}_{1}})}}\left\lbrack {\log_{2}\frac{\rho}{t}{{\mu\lambda}_{i}^{2}\left( {S_{2}{HS}_{1}} \right)}} \right\rbrack^{+}}}}},$ where μ depends on (S₁,S₂) and satisfies

${\sum\limits_{i = 1}^{\min{({L_{t},L_{r}})}}\left\lbrack {\mu - \frac{t}{\rho\;{\lambda_{i}^{2}\left( {S_{2}{HS}_{1}} \right)}}} \right\rbrack^{+}} = {t.}$

FFT-based Method (FFTS)

Similar to the FFTS in diversity transmission case, the FFT matrix is inserted in the RF chains together with the selection switch. The optimum selection with FFT can deliver a capacity of

$\begin{matrix} {{C_{FFTS} = {\max\limits_{S_{1} \in S_{t}}\;{\max\limits_{S_{2} \in S_{r}}{\sum\limits_{i = 1}^{\min{({L_{t},L_{r}})}}\left\lbrack {\log_{2}\frac{\rho}{t}\mu\;{\lambda_{i}^{2}\left( {S_{2}\Phi_{r}H\;\Phi_{t}S_{1}} \right)}} \right\rbrack^{+}}}}},} & (16) \end{matrix}$ with the constraint

${\sum\limits_{i = 1}^{\min{({L_{t},L_{r}})}}\left\lbrack {\mu - \frac{t}{\rho\;{\lambda_{i}^{2}\left( {S_{2}\Phi_{r}H\;\Phi_{t}S_{1}} \right)}}} \right\rbrack^{+}} = {t.}$

Phase-Shift & Selection Method (PSS)

For PSS method, we use the same notations as before

$C_{PSS} = {\max\limits_{F_{1} \in F_{t}}\;{\max\limits_{F_{2} \in F_{r}}{\sum\limits_{i = 1}^{L}\left\lbrack {\log_{2}\frac{\rho}{t}\mu\;{\lambda_{i}^{2}\left( {F_{2}{HF}_{1}} \right)}} \right\rbrack^{+}}}}$

Partial CSI

In the case that the channel is only partially estimated, the invention provides an alternative method. In this case, there is no knowledge about the instantaneous channel state, but rather only about the state of the channel averaged over a time interval, e.g., typically, but not restricted to, several milliseconds to several hundreds of milliseconds. This case typically occurs for the transmitter in a frequency-division duplexing system,

In that case, the goal of the transformation is to transform the transmit signal in such a way that on average the signals after the selection process has good properties, e.g., high capacity.

The expected value capacity that can be achieved that way is

${E\left\{ C_{PSS} \right\}} = {\max\limits_{F_{1} \in F_{t}}\;{\max\limits_{F_{2} \in F_{r}}\;{E{\left\{ {\sum\limits_{i = 1}^{L}\left\lbrack {\log_{2}\frac{\rho}{t}\mu\;{\lambda_{i}^{2}\left( {F_{2}{HF}_{1}} \right)}} \right\rbrack^{+}} \right\}.}}}}$ The expected value of the capacity is only one possible criterion. Other criteria result in different values for the capacity matrix.

Effect of the Invention

FIG. 2 compares the SNR gain of the four methods described above as a function of antenna spacing with diversity transmission. The parameters in the system are t=8, r=8, L_(t)=1, L_(r)=2. For large correlation, i.e., small antenna spacing, the FFTS method in the invention has a considerably larger SNR gain over the traditional HS-MRT/MRC method. In general, regardless of the correlation level in the channel, the PPS method in the invention can outperform both prior art selection methods. In fact with two or more RF branches, the PSS method can achieve the SNR gain of a full-complexity method.

FIG. 3 displays the 5% outage capacity of the four methods for multiplexing transmission. The system parameters are t=3, r=8, L_(t)=3, L_(r)=3, p=20 db. A similar conclusion can be drawn from the plots: the FFTS shows a performance improvement for strongly correlated channels while the PPS works well for both strong and weak correlations.

This invention is described using specific terms and examples. It is to be understood that various other adaptations and modifications may be made within the spirit and scope of the invention. Therefore, it is the object of the appended claims to cover all such variations and modifications as come within the true spirit and scope of the invention. 

1. A method for processing radio frequency (RF) signals in a multi-antenna system, comprising: generating L_(t) input data streams in a transmitter with t transmit antennas, where t is an integer, and where L_(t) is an integer less than or equal to t; modulating the L_(t) input data streams to RF signals; switching the RF signals to t RF branches; applying a phase-shift transformation to the RF branches by a t×t matrix multiplication operator Φ₁, whose output are t RF signals, and in which entries of the matrix Φ₁ have constant modulus phase-only terms and the modulus phase-only terms adapt to an estimate of an average channel state; transmitting the t RF signals over a channel by the t transmit antennas; receiving the transmitted signals in a receiver with r receive antennas, where r is an integer; applying a phase-shift transformation to the r RF signals by a r×r matrix multiplication operator Φ₂ to generate r streams; selecting L_(r) signal streams from the r streams, where L_(r) is an integer less than or equal to r; demodulating the L_(r) signal streams; and processing the demodulated L_(r) signal streams in baseband to recover output data streams corresponding to the input data streams.
 2. The method of claim , in which each of the L_(t) input data stream has a weight, and further comprising: summing the L_(r) weighted data streams before the demodulating and decoding.
 3. The method of claim 1, in which the L_(t) input data streams are generated by a space-time block coder.
 4. The method of claim 1, in which the L_(t) input data streams are generated by a space-time trellis coder.
 5. The method of claim 1, in which the input data streams are space-time layered structures.
 6. The method of claim 1, in which t=L_(t), and the matrix Φ₁ is an identity matrix.
 7. The method of claim 1, in which r=L_(r) and the matrix Φ₂ is an identity matrix.
 8. The method of claim 1, in which the phase-only terms adapt to an estimate of an instantaneous channel state.
 9. The method of claim 1, in which entries of the matrix Φ₂ have constant modulus phase-only terms.
 10. The method of claim 1, in which entries of the matrices Φ₁ and Φ₂ have constant modulus phase-only terms.
 11. The method of claim 1, in which the matrix Φ₁ is a fast Fourier transform matrix.
 12. The method of claim 1, in which the matrix Φ₂ is a fast Fourier transform matrix.
 13. The method of claim 1, in which the matrices Φ₁ and Φ₂ are fast Fourier transform matrices. 